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SUMMARY 


As aircraft Income larger and lighter due to design requirements for 
increased payload and improved fuel efficiency, they may also become much 
more flexible. For highly flexible vehicles, the handling qualities may not be 
accurately predicted by conventional methods. This study applies two analysis 
methods to a family of flexible aircraft in order to investigate how and when 
structural (especially dynamic aeroelastic) effects affect the dynamic 
characteristics of aircraft. The first type of analysis is an open-loop modal 
analysis technique. This method considers the efTect of modal residue 
magnitudes on determining vehicle handling qualities. The second method is :• 
pilot- in-the-loop analysis procedure that considers several closed-loop system 
characteristics. Both analyses indicated how dynamic aeroelastic effects can 
cause a degradation in vehicle tracking performance, based on the evaluation of 
some simulation results. 

This report is divided into two volumes. Volume I consists of the 
development and application of the two analysis methods described above. 
Volume II consists of the presentation of the state variable models of the 
flexible aircraft configurations used in the analysis applications, mode shape 
plots for the structural modes, numerical results from the modal analysis, 
frequency response plots from the pilot-in-thc-loop analysis and a listing of the 
modal analysis computer program. 


* 
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CHAPTER I 
INTRODUCTION 


Usually, the attitude dynamics and handling qualities of aircraft are 
defined in terms of rigid-body modal characteristics. For example, the 
frequency and damping of the short-period and phugoid modes are used for 
handling qualities specifications of the longitudinal dynamics of aircraft [1]. 

This is possible, not because the aircraft are actually rigid, but because they 
are “rigid enough” so that structural effects can be ignored. It has been shown, 
however, that this approach may not be very accurate for aircraft with 
significant amounts of structural flexibility [2]. Since, in the future, aircraft 
will become larger and lighter due to design requirements of increased payload 
and improved fuel efficiency, they may also become much more flexible. In 
addition to the rigid-body modes, flexible aircraft have aeroelastic-structural 
modes which may significantly affect their dynamic characteristics. Analysis of 
the dynamics of these aircraft, without considering the contribution of the 
structural modes, would be inaccurate. Any use of such an analysis, in flight 
control designs for example, could produce poor, if not disastrous results [3]. 

At present there is no universally accepted way to predict the handling 
characteristics of an aircraft in which structural flexibility is significant. 

Further, there is a need to describe qualitatively, the significance of structural 
effects. The goal of this research, then, is to address the questions of when and 
how do structural effects (especially dynamic aeroelastic effects) significantly 
affect the dynamic characteristics of aircraft? Answering these questions is the 
first step in developing a systematic approach to analyzing flexible aircraft 
handling qualities and synthesizing appropriate flight control laws. 

This report is divided into the following chapters that present the 
development and application of the analysis “tools”. Chapter 2 uses pole-zero 
plots and transfer functions of flexible aircraft to provide background on why 
structural effects can be significant and therefore explains why they need 
consideration. Chapter 3 presents the family of vehicle configurations which 
will be used throughout the analysis. In Chapter 4, an open-loop analysis 
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CHAPTER U 
BACKGROUND 


When the vibrational frequencies of an aircraft structure are large 
compared to the frequencies of the rigid-body modes, the effect of the flexible 
modes on the overall dynamic response of the aircraft is small. This is the 
situation for most aircraft and, as will be seen, allows the dynamics to be 
accurately modeled by the rigid-body modes only. However, as the frequencies 
of the structural modes become lower, the effect of these modes on the 
dynamics can become significant. 


For example, consider the attitude response of an hypothetical aircraft due 


to elevator deflection by studying the transfer function 


JfeL 


«•) 


, or the pole-zero 


plot corresponding to this transfer function. Figure 2.1 shows a typical pole- 
zero plot of the longitudinal attitude response transfer function where the poles 
and zeros of the phugoid and short-period modes and the first few structural 
modes are included. Four poles and two zero may be considered to be 
associated with the phugoid and short-period modes. Typically, the poles are 
complex conjugates and the zeros are real. Note also that there is a pole-zero 
“dipole" associated with each of the struct^al modes. The poles and the zeros 
are complex for these modes. 


Although a real aircraft, like any structure, has an infinity of vibrational 
modes, for ease of discussion the example used here will consider only one of 
the structural modes. To further simplify the discussion, the phugoid mode 
will also be omitted, that is, the “short-period approximation" is invoked. 

The pole-zero plot simplifies to Figure 2.2 when the above simplifications 
are applied. The transfer function associated with this simplified case appears 
in Equation (2.1). 
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Figure 2.1 

Pole- Zero Plot of Typical Flexible Aircraft 





Figure IS 

Poie-Zero Plot of Simple Example 
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fl(s) _ (s-z sp )(s-z 1 )(s-Z|) 

M s ) S (s-p S p) (s-p sp ) (s-p,) (s-p,) 

- where (“) denotes the conjugate of ( ). Equivalently, the transfer 
function for pitch rate due to elevator input is, 


fl(s) _ (s-z, p )(s-zi ) (s-z l ) t 

W (s-p sp )(s-p gp )(s-p,)(s-p,) 

The pitch-attitude- rate response of the aircraft due to an impulsive input 
is therefore, 


0(s) = 


(s-z 8p )(s-Zt)(s-z,) 
(s-p*p) (s-p sp ) (s— p,) (s-p,) 


(2.3) 


The following form of the attitude-rate response results from the partial 
fraction expansion of Equation (2.3) and transformation into the time domain. 


0(t) = R sp e p * 1 + R, p e p ^ + R^* 1 + R,e p,t (2.4) 

Here Rj is the residue associated with pole pj, and R; is its conjugate. For 
complex poles, the residues ire complex numbers with magnitudes that 
determine the degree to which each mode contributes to the overall response. 
Therefore, the significance of an individual mode in the dynamics of the 
aircraft is represented by the residue magnitude of that mode. This can be 
illustrated by writing Equation (2.4) in the form, 

0(t) = 2 1 R gp | e ff,pt cos(w ep t + * sp ) 

+ 2 1 R,| e ff,t co$(u>,t + $,). (2.5) 


i 


- where Pi = <J\ + jw; , 
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4> ; = tan 1 


Im(R i ) 


R; I = 


Re(Rj) 

[Re(Rj)] 2 + |Im(Ri)] 2 


m 


The residue magnitudes can be interpreted geometrically in terms of the 
poles and zeros of Figure 2 2 by considering the following relations [4], 


and, 


p I — ^gp*~Pgp I z rP sP ' P$p 

sp | Psp-P*p|i PrP.pM PrP.p| 

R | - I 2 i~PiM z i~PiM VPiI 
R sp I — — 

Psp Pi ’ Pgp Pi ‘ Pi Pi 


(26) 


(2.7) 


and from symmetry of the pole-zero plot, j R-, | = | Rj| . Here | '| denotes the 
magnitude of a complex number and so | z a — pj| is the distance from the point 
z, to the point pj. 

If the hypothetical aircraft was fairly rigid, the poles and zeros associated 
with the structural mode would be far from the origin when compared to the 
poles and zero associated with the short-period mode. This implies that, 


and, 


Zl-Pgp 1 

“ 1 PrPgpl 

(2-8) 

? i~p*p 1 

“ I PrPgp I • 

(2.0) 


In this case, the expression for the residue magnitudes associated with the 
short- period mode can be simplified by the effective cancellation of terms 
involving z l and p, from the numerator and the denominator so that the 
short-period residue magnitudes are relatively independent of the structural 
mode, or, 
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R r- 

ll, sp — 


Z 8P Psp 


Psp Psp 


( 2 . 10 ) 


Similarly, for a fairly rigid aircraft, the zero near a pole associated with the 
structural mode is much closer to that pole than any other pole or zero of the 
system. Also, the distance from that pole to its complex conjugate is 
approximately equal to the distance to the zero associated with that conjugate 


pole (i.e. 

the zero of the pole-zero dipole). 

These two statements imply that, | 

* 1 


I VPi | 

« 1 Psp-Pl 1 . 

(2.11) j 


1 Z|"Pl | 

« 1 i> S p“Pi 1 » 

(2-12) j 

and, 



f 

3 


i vpi 1 

- | pr Pi | • 

(2.13) ! 


Here, (he residue magnitudes associated with the structural mode can be 
simplified, using Equation (2.13), to obtain the following expression. 


, ,V | .,l l ! r >, l <0 . 


P*p-Pi | | I VP I 


(2.14) 


It is clear, by applying Equations (2.11) and (2.12), that the structural mode 
residue magnitude is much less than unity. Since the short-period mode residue 
magnit ude is on the order of unity, it is obvious that 


| -l|| « | ■ <2.151 

As a result of this discussion, two conclusions can be drawn concerning 
fairly rigid aircraft. First, the structural modes have little affect on the degree 
to which the rigid-body modes contribute to the dynamic response of the 
vehicle. And second, the contribution of the structural modes to the dynamic 
response of the vehicle is insignificant compared to the affect of the rigid-body 
modes. Therefore, the longitudinal attitude-rate impulse response of the 
aircraft can be accurately approximated by the following expression. 
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«t)»2| R gp | e^cosfwgpt + * sp ) (2.16) 

The implication of this expression is that, when it is valid, the dynamics of the 
aircraft are determined almost entirely by the values of the rigid-body poles 
and zeros. 

However, if the amount of flexibility is significant enough so that 
Equations (2.8), (2.9), (2.11) and (2.12) are no longer valid, the conclusions 
above will no longer apply. In this case, the degree to which the rigid-body 
modes contribute to the overall response will depend somewhat on the 
characteristics of the structural modes. In addition, the contribution of the 
structural modes to the response may be significant. 


CHAPTER HI 

EXPERIMENTAL DATA BASE 


A set of aircraft dynamic models, one of which is similar to the B-l 
bomber, were available from a previous study [2]. The B-l is a large aircraft 
with a reasonable amount of structural flexibility. Figure 3.1 is a sketch that 
depicts the geometry of the aircraft that corresponds to all the vehicle models 
to be considered. The models represent a family of aircraft similar to the B-l 
that differ only in their amount of structural rigidity, quantified in terms of the 
invacuo-structural vibration frequencies. The configurations can be described 
physically as vehicles with identical geometries but made of different materials 
so that the vibration frequencies are changed while the vibration mode shapes 
remain unchanged. 

The mathematical models of the aircraft include two structural modes 
which correspond to the first fuselage bending mode and the second fuselage 
bending mode. The mode shapes which correspond to these acroelastic- 
structural modes can be found in Appendix A.4. The family of configurations 
were generated by parametrically varying the invacuo-structural frequencies of 
the two structural modes. Table 3.1 summarizes the eight configurations, 
listing their eigenvalues and the invacuo-vibration frequencies of the two 
structural modes. 

Notice that for configurations 6, 7 and 8 the second aeroclastic mode is 
slightly unstable due to negative aerodynamic damping. The original 
simulation study [2] involved considering the effect of neutrally stable modes on 
vehicle dynamics. To study this effect, the very slightly unstable (i.e. 
effectively neutrally stable) configurations 6, 7 and 8 were developed. 

The complete mathematical model of the eight configurations in state 
variable form corresponding to Equation (3.1) can be found in Appendix A.3. 
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iv = Av* v + B v ii + D v w 


(3.1) 


- where x v is the vector of vehicle states and Ay, By and D v are system 
matrices, u is the vector of control inputs and & is the vector of disturbances. 


Table 3.1 

Summary of Data Base Configurations 



!NV> 
MODE FI 

\CUO 

?EQ’S (H*P 

AE 

M 

ROELASTIC VEHICLE 
IODE EIGENVALUES 

CONFIG 

MODE 

1 

MODE 

2 

PHUGOID 

SHORT 

PERIOD 

MODE 

1 

MODE 

2 

1 

2.18 

3.37 

•0.0015 

ij0.067 

-1.5 

±j2.37 

-0.66 

±jl3.3 

-0.46 

±j21.3 

2 

1.46 

3.37 

0.001 

dbj0.053 

-1.35 

±j2.2 

-0.73 

±j8.76 

•0.46 

±j21.3 

3 

0.97 

mm 

-0.08 ; 
0.095 

-0.9 
dr j 1 .5 

-1.11 

dtj5.7 

•0.46 

±j21.3 

4 

2.18 

0.76 

-0 .13 ; 
0.15 

-1.06 
L ±1L1 

-0.7 

±jl3.3 

^>.53 

±j5.9 

5 

1.86 

I 86 

•0.001 

±j0.049 

-1.4 

±j2.17 

•0.86 

±jll.7 

-0.12 

±jll.6 

6 

1.10 

1.10 

-0 .15 ; 
0.18 

-0.95 

ijO.97 

•1.31 
±j7 16 

0.057 

±j7.0 

7 

1.63 

1.55 

0.001 

±j0.017 

-1.32 

±j2.0 

-1.1 

±jl0.2 

d-H 

8 

1.70 

1.48 

0.0013 

±j0.012 

-1.3 

±]2.0 

•1.08 

d:jl0.3 

0.085 

L-±M ■ 


In the previous study [2], the above vehicle configurations were used in a 
fixed based, laboratory simulation involving longitudinal tracking of a low 
frequency command signal. A cathode ray tube was used to display the 
following variables; commanded attitude angle, and vehicle attitude angle, 
measured at the cockpit location. The vehicle attitude angle is the pitch 
attitude measured, for example, by a gyro located at the cockpit and differs 
from the rigid vehicle pitch Angle, 0 R , by the contribution of the local 
structural deflections. This effect is illustrated in Figure 3.2 and defined by 
Equation (3.2). 









(3.2) 


«l(t) = w<)- i>' (I, ll,(t| 

i=l 

- where l x is the cockpit location measured from the center of gravity, <p' 
is the mode slope of the ith elastic mode and is the generalized coordinate of 
the ith elastic mode [5j. 

The above information was displayed to the pilot by meane of a visual 
display similar to the one depicted in Figure 3.3. His task was simply to 
minimize the error between the commanded and the indicated attitudes angles. 

Three types of data were collected from the simulations : 1) rms tracking 
error (taken over a 120 second run for each case); 2) Cooper-Harper (6] pilot 
rating in the tracking task *; and 3) pilot comments. A summary of these 
results can be found in Table 3.2. These results indicate hat flexible aeroelastic 
effects significantly affected pitch attitude tracking performance. 


Table 3.2 

Summary of Tracking Error, Pilot Rating and Pilot Comment* 


CONFIG 

RMS Tracking 
Error (d^g) 

RMS Pilot 
Rating 

Pilot Comments 

1 

1.13 

1.6 

Good; no problem 

2 

1.05 

2.0 

little oscillation; slight 
control response lag 

3 

5.67 

5.9 

difficult; PIO problem; 
extreme response lag 

4 

1.90 

3.1 

little more difficult than Cl; 
sluggish attitude response 

5 

1.51 

2.0 

pretty good: same as C2 

e 

7.57 

6.7 

severe oscillation; vir- 
tually uncontrollable 

7 

1.48 

2.3 

not difficult; annoying 
oscillation 

8 

1.16 

1.9 

not difficult; little 
oscillation but could ignore 
it and By rigid body 


note: These reaulta are for 4 pilots with 2 runs per pilot. 



i 
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I 

Before introducing the analysis methods, *he simulation results, which are 
summarized in Figures 3.4 and 3.5, will be reviewed. It is clear that aeroelastic 


f The Cooper-Harper rating stale is a subjective rating (1 being best and 10 worst) used 
to describe vehicle handling qualities in various tasks. (G| 
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Figure 3.3 

Simulation Visual Display 
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Figure 3.4 

Simulation Tracking Error* 


pilot rating (cooper-harper) 


level 1 
level 2 
level 3 

- good 

- fair 

- poor 


«• 

I 

0 

level 3 

Q 

1 

level 2 


1 2 3 4 5 6 7 8 

configuration 


Figure 3.S 

Simulation Pilot Ratings 




effects significantly affected vehicle dynamics. By merely varying invacuo- 
structural frequencies, the dynamics changed so drastically that two 
configurations (3 and 6) received Level 3 ratings while the others received Level 
I ratings. Once again the question to be answered is, “when and how do these 
aeroelastic effects affect aircraft dynamics?” 

Note that the “rigid-body” phugoid and short-period eigenvalues alone 
give little insight into the effect of reducing the invacuo-structural frequencies 
(see Table 3.1). For example, Configuration 3 has a much higher (worse) 
Cooper-Harper rating and larger tracking errors than Configuration 4 despite 
the fact that Configuration 3 has a more stable phugoid mode and only a 
slightly higher frequency associated with the short-period mode. In addition, 
Configurations 3 and 4 have similar lowest- frequency aeroelastic mode 
eigenvalues. Based on this, one might predict that Configuration 4 should be 
worse than Configuration 3, which is contrary to the simulation results. Thus, 
the eigenvalues alone do not completely capture the actual dynamics of a vehicle. 
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CHAPTER IV 

OPEN-LOOP MODAL ANALYSIS 


A modal analysis of the family of aircraft presented in Cnapter 3 will be 
presented. The results of this analysis will be used to help explain some of the 
findings which were obtained in the simulation of those dynamic configurations. 
These results will also be used to attempt to answer the questions posed earlier 
• namely, when and how do structural effects significantly affect the dynamics 
of aircraft? 


Modal Analysis 

Consider the vehicle modeled in the state variable form, 
i = As. + Bu + Dw 


(4.1) 


y = Cx. + Eh + Fw 

Here \ is a vector of vehicle states, y is a vector of outputs and u and s y arc 
vectors of control inputs and disturbances, respectively. 

One may diagonalize the system using the modal transformation (4,5), 
where the modal matrix, T, is formed from the eigenvectors of A ( assuming 
distinct eigenvalues ), so that 


T $ lni.iis kJ 


(4.2) 



Here % is the eigenvector associated with the ith eigenvalue of A. In terms of 
the modal states, the system dynamics are, 

g = Ag + Bg + Dw , 


(4.3) 


y = Cg + Eg + Fw . 

Here g is the vector of modal coordinates, A ^ T~'AT (diagonal), B ^ T _1 B, 
D T~*D and C = CT. The matrices B, C and D are called the modal 
controllability, disturbability and observability matrices, respectively. With 
proper vehicle state definitions, elements of these matrices indicate how 
controllable, disturbable and observable each mode is, with respect to the 
inputs and outputs. Each element of the modal controllability matrix, for 
example, is a relative measure of how much the associated control input 
contributes to the response of the corresponding mode. If the magnitude of one 
element of the modal controllability matrix is small compared to the 
magnitudes of the other elements of the controllability matrix then the mode in 
question is “relatively uncontrollable” from that input. 

At this point, to simplify the development, one may redefine the input 
vector as, 

a £ [jj]. H.4) 

Substituting this into Equation (4.3) results in, 
g = Ag + Bu 


(4.5) 


X = Cg + ES, 


where Q ^ 


B 

D 


and £ 




The diagonal property of A is especially useful when considering the 
system in the frequency domain, or, 
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[sI-A]gH = Bu(s), 


jl(s) - Cg(s) + Eu(s). 


(4.6) 


Since the identity matrix, I, and A are diagonal and square, the inverse of 
[si- A] is diagonal as well. By multiplying the first of Equations (4.6) by 
(sI-A) 1 and substituting into the second of Equations (4.6), the following 
matrix equation for the outputs y results. 


y(s) =C[sI-A] 'Bli(s) + Eu(s) 


(4.7) 


My writing the elements of C as Ojj (i the row and j the column), the 
elements of B and E as bj: and e^, respectively, and by using the fact that 
[ si — A 1 1 is diagonal, the transfer function for the ith output due to the jth 
v,(s) ' 


input . 


"j(^) 


, can be written as, 


y,< s > 

'lj(s) 


= V 
k = l 


f .k • hi 

S X L 


(4.8) 


Here n is the number of system states and X k is the kth eigenvalue of the 
system. 

From Equation (1.8) the impulse response of yj can be obtained by 
assuming itj to be an impulse and taking the inverse Laplace transform. The 
impulse response of y-j becomes, 

>i(t) = U c ik ^kj «*X|»(X k t ) + e ijt (1.0) 

k = l 


Note that the values of c^ • by are the values of the residues associated with 
mode k (k = l,2,...,n) for the ith output (y 4 ), when the system is excited by the 
jth input (Uj). This relationship results from the definition of a “residue" [4] 
and Equations (4.8) and (4.0) so that, 
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\ 


R 


k (y, u,) c ik '^kj- 


(4.10) 


Therefore, the information reflected in the controllability, disturbability and 
observability of a mode is also completely contained in the modal residues. 
Equation (4.9) can therefore be written as, 


yi(t) = E R k ex P( x k l ) + e ij- (4-11) 

k = l 

- where n is the number of system poles. By representing the residue in 
terms of its magnitude and phase and combining terms involving complex 
conjugates, Equation (4.11) can be written as, 

yi(t) = E 2 | R k| e'* kt cos(w k t + * k ) + e^. (4.12) 

k = l 

- where m is the number of modes of the system (i.e. m = n/2) and all 
eigenvalues are assumed to be complex for ease of discussion. As discussed in 
Chapter 2 and clearly from the above equation, the magnitude of the residue of 
a mode is a direct measure of the contribution of that mode to the dynamic 
response of the vehicle. From Equation (4.12), the relative importance of each 
mode to a given response can be determined by inspection. By numerically 
implementing this analysis, it is possible to investigate how higher order modes 
directly affect the dynamic response of an aircraft by comparing impulse 
residue magnitudes. 

The results of the modal analysis as developed above include the 
eigenvalues and eigenvectors of the system. The eigenvectors are used to 
determine which mode of the dynamics is associated with each generalized 
coordinate of the system. It should be noted that when the term “rigid-body 
mode” is used it means the system mode whose eigenvector reflects significant 
participation of the rigid-body states (e.g. attitude, attitude rate and angle of 
attack). Similarly, the term “elastic mode" is used to mean a system mode 
whose eigenvector reflects significant participation of the elastic states (i.e. ifc, 
»)i). In the application to flexible aircraft, there are no truly rigid-body modes 
or purely structural modes due to aeroelastic coupling. 

The manner in which the modes of the dynamics associated with the 
generalized coordinates are identified can best be explained through example. 
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Consider the data presented in Table 4.1 which represents the two rigid-body 
modes of the Navion aircraft [7j. The magnitude of the element associated 
with attitude rate, 0, in the first eigenvector is larger than the elements 
associated with the other states. This indicates that the mode associated with 
the first eigenvector primarily contains attitude rate and therefore, corresponds 
to the short-period mode. Similarly, for the second eigenvector, the magnitude 
of the element associated with forward velocity perturbation, u, is larger than 
the other elements and so, that mode corresponds to the phugoid mode. This 
technique was used to identify the vehicle modes for the configurations used in 
this study. 


Table 4.1 

Summary of Navion Longitudinal Dynamics 


Flight ConditioD : Sea Level 

U = 176.0 ft/sec 
W ■ 1.84 ft/see 

State Vector : ^ 1 u, w, 0, 0 1 


units 

eigenvalues 

m/mmm 

■Mol 

1 

Bill 

sec -1 

eigenvectors 


1.0 

ft-sec“* 



ft -sec" * 

1.0 

0.143 
± jO.009 

10 -»M 

•0.193 
± jO. 190 


10" 3 rad 


Notice that, in the example, the units of the elements of the eigenvector 
are not the same. Two of the elements have units of feet per second and the 
other two have units of 10" 3 radians per second and 10 -3 radians, respectively. 
This set of units enable the modes of the system to be readily identified. Since 
the magnitude of an element of an eigenvector is dependent on the units 
selected for the system states, proper choice of the units can aid in identifying 
the modes of the system. In general, the states of a system (and therefore the 
elements of the eigenvectors) do not have the same physical units. The units 
can be changed, however, by applying a similarity transformation to the state 
variable representation of the system. It can be easily shown that a similarity 
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transformation does not affect the eigenvalues or residues of a system. As a 
result, a similarity transformation can be applied to the system without 
altering the modal analysis results. Appendix A. I presents the development of 
the transformation used in this study to change the units of the states and 
eigenvectors to aid in identifying the system modes. 

In addition to the eigenvalues and eigenvectors, the modal controllability, 
disturbability and observability matrices are available from the modal analysis. 
With proper selection of units (using a similarity transformation), these can be 
used to gain further information concerning dynamic relationships between the 
control and disturbance inputs and the system modes, and between the outputs 
and the system modes as described previously. 

The last and most useful of the analysis results are the modal impulse 
residues, Rj. The impulse residues are useful since, as noted previously, they 
are a direct combination of the observability and controllability (or 
disturbability) of a particular mode. The magnitude of the impulse residues 
will be used extensively in this study. 


Vehicle Model 

The vehicle models considered were those used by Yen [2]. They consist of 
linear equations of the form • 


*R 

*E 


Ar 

A c ' 


Ac 

A e 


*R 

*E 


+ 




(4.13) 


The vehicle states include some states which are rigid-body degrees of freedom, 
X R , and others which are the structural degrees of freedom, * E . Together these 
form the partitioned state vector. The system matrices are also partitioned to 
be consistent with the state vector partition. The sub-matrices with the 
subscript R are those associated with the rigid states and the sub-matrices with 
the subscript E are those associated with the elastic states. The sub-matrices 
Ac and relate the cross-coupling between the rigid states and the elastic 
states. 

The modal analysis procedure could be accomplished with the Oexible 
vehicle model described in Equation (4.13) but the results would not be the 


most meaningful. This is due to the fact that interpretation of the modal 
analysis requires extensive use of characteristic parameters that influence the 
impulse response of the system (i.e. modal controllability, modal observability, 
residues, etc.) and since an impulse input is unrealistic, these characteristic 
parameters are unrealistic. Since an impulse is physically unrealizable, impulse 
responses of an aircraft are unrealistic and do not reflect the dynamics of a 
vehicle in actual flight. In order to obtain meaningful results from the analysis, 
inputs should represent important aspects of the actual pilot commands and 
atmospheric turbulence. 

An impulse has infinite bandwidth and cannot be produced by any 
physical system. A pilot may try to produce an impulse but, due to his limited 
bandwidth, cannot achieve it. What results is a type of "realistic pulse” input 
that can be approximated by treating the pilot as a low pass filter (i e. a first 
order lag). The input to the filter is an impulse and the resulting output is a 
finite bandwidth pulse which approximates what a pilot is capable of 
producing. Equation (4.14) is the state space representation of a low pass filter. 


^ p' A p x p ^*p9p 


(4.14) 


The scalar x p is the "realistic pulse”, t| p is the impulsive input, Ap = — — , 

l 

where r p is the time constant of the filter, and G p is — so that the Bode gain 

r p 

of the filter is unity. By using x p as the input to the vehicle model and an 
impulse as the input to the filter, a realistic response is obtained for the 
pilot/vehicle system. 


The importance of using the low pass filter to obtain meaningful results 
can be clearly shown by a simple example. Consider a system with two states 
in modal coordinates - 


i = 


-l 

o 



l 

l 


u. 


( 4 . 15 ) 
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The response of this system to an impulse input, u = 6(t), is, 
y(t)| u =<(t) = e 1 + e ,00t • 


(4.16) 


Now consider the same system and let the input, u, be represented as, 

u = -lOu + 10 f| . (4.17) 

That is, u is the output of a low pass filter with a time constant, r p = 0.10 sec*. 
The system can be reorganized as, 




-1 

0 

1 


0 

i 

= 

0 

-100 

1 

w + 

0 

u 


0 

0 

-10 

luj 

10 


(4.18) 



L 


The response of this system to an impulse input, »j = 6(t), is, 

=(l.ll)e' t 

+ ( 1.00) e* 101 + (0.11)e' ,00t . (4.10) 


Note the difference between the two responses, Equations (4.16) and (4.10). The 
unrealizable response (Equation (4.16)) indicates that both modes contribute 
equally to the '"rerall response, in terms of residue magnitude. The realizable, 
filtered response (Equation (4.10)) indicates that the fast system mode 
(X =-100) has a much less significant contribution to the overall response than 
does the other original system mode (X = - 1). The obvious conclusion is that 
an impulse input to the system excites the fast mode, but cannot be excited as 
much by the limited bandwidth filtered impulse. Therefore, the modal analysis 


t A time const sat of 0.10 seconds is consistent with human bandwidth limitations. 
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should be performed on the system which includes the low pass filter that more 
accurately reflects the true inputs that are expected. If the filter is not used, 
the modal analysis may indicate that certain high frequency modes significantly 
contribute to the vehicle response when they actually may have insignificant 
effects. 

An argument similar to that used for the pilot commands can be used to 
justify describing the disturbances produced by the atmosphere in a similar 
way. Since atmospheric turbulence is an important aircraft disturbance and 
turbulence is random in nature, ihcse disturbances are modeled stochastically. 
One commonly used disturbance model is the Dryden Gust Model (7). The 
Dry den Model may be expressed in matrix form as, 


— l A g,' n, gl » 

- A g ij + G g »; g . 


(4.20) 


The gust state n g is the angle of attack induced by a vertical gust and a fi| is an 
additional gust state which is necessary to obtain the proper frequency 
character of the gust model. The system matrices A g and G g (given later) 
provide the proper characteristics of the random gust response when the 
“white" noise, tj t , is the input. 

Hy combining the “pilot equation" and the Dryden Gust Model equation, 
(Equations (4.14) and (1.20)), with the aircraft, (Equation (4.3)), an augmented 
flexible aircraft mathematical model is formed. The resulting flexible aircraft 
model is then composed of a combination of the system matrices from the pilot 
equation, the gust equation and the flexible vehicle equation, (Equation (4.22)). 




(4.21) 



> 

O 

O 


Gp 


0 

i = 

0 A ( 0 

By Dy Ay 

•* + 

0 

0 

*p + 

G* 

0 


(4.22) 
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Finally, the proper choice of aircraft responses, y, is critical for obt: ning 
meaningful results from the analysis. If the wrong outputs are selected, 
erroneous conclusions may be drawn. This point is emphasized so that the 
reader b aware that a great deal of engineering judgement must be used in 
choosing the proper outputs. Understanding the physics of the problem is 
necessary to obtain meaningful results. 

Once the outputs of interest are chosen, linear matrix output equations are 
formed so that, 

y = Ci + E»7 p + F»/ g (4.23) 

- where y is a vector of outputs. The modal analysis may now be performed, 
using Equations (4.21), (4.22) and (4.23) as the complete system. 


Application To Data Base Configurations 

The modal analysis method was implemented in a computer program. A 
listing of this program appears in Appendix A.6. As a result of the modal 
analysis, several quantities of interest are readily available. In addition to the 
modal impulse residues, the modal eigenvalues, the modal eigenvectors, and the 
modal controllability, dbiurbability and observability matrices are all easily 
obtained. 

The "chicle models used in Yen's simulation were extended to include two 
additional structural modes. The additional modes were t.«e second and fourth 
lowest frequency modes of the baseline vehicle (i.e. Configuration 1), thus 
increasing the model to include the four lowest frequency structural modes. 

The mode shapes of the additional modes indicate that they are primarily 
symmetric wing bending modes. The shapes of these modes can be seem in 
Appendix * 4. These modes could be important in the gust responses and will 
be considered later in the analysis. 

The total state vector includes the standard rigid-body degree of freedom 
(i.e. perturbed forward velocity, u ; angle of attack, o; rigid-body pitch attitude 
and attitude rate, 0 R and 0 r) and the generalized coordinates of the four 
structural modes, (i. the generalized deflections, ^ ; and the generalized rates 
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ij,). The total state vector is defined as, 
x v T ^ Ur | x e t ] 


= [or, 6? r , U, <? R | tit | f) 4 l. 


(4.24) 


The system matrices of the vehicle configurations using this state vector can be 
found in Appendix A.3. 

The output parameters were chosen to include rigid-body flight path angle 
('/), (Equation (4.25)). total-elastic (0 T ) and rigid-body (0 R ) pitch attitude 
angles, (Equation (4.26), Figure 4.1), total-clastic (fy.) and rigid-body (0 R ) pitch 
attitude rates, (Equation (4.27)), and normal acceleration at the cockpit (n,), 
(Equation (4.28)). 


7r ~ °R~ q r > (rad) 

(4.25) 

h ~ °R~)2 1 lM) < i > ' < (rad) 
1=1 

(4.26) 

*T = «R-i;iMlVilU .(rad/sec) 
i = l 

(4.27) 

". = 7 ieoi + iA-i«m(U! Os's) 

? i=l 

(4.28) 

-where g = gravitational acceleration, (ft/sec 2 ) 


U 0 = cruise velocity, 949 (ft/sec) 


l x = distance between c.g. and cockpit, (ft) 




\tange nt to 
'body at 

cockpit location (1 ) 


Figure 4.1 

Rigid and Elastic Pitch Angles 
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The above parameters constitute what was judged to be the significant 
responses in longitudinal attitude dynamics. Total-elastic and rigid-body pitch 
attitude angle and pitch rate are used extensively by the pilot to control the 
vehicle and evaluate its performance. In fact, in the simulation study, the 
pilot’s task was to minimize the error between a commanded attitude, 0 C , and 
the vehicle attitude, 0 T . This implies that 0? and as well as 0 r and 0R are 
of extreme importance in pitch attitude tracking. Normal (or plunge) 
acceleration is another significant response of the vehicle from the aspect of 
ride quality, but of course was not a factor in the fixed-base simulation. 

Note that the equation for n, , (Equation (4.28)), is not an explicit 
function of the states in Equation (4.24), but is a function of the state 
derivatives. It is, therefore, an implicit function of the system states and 
control deflections. By using the state equations (Equation 4.1), n, can be 
written as an explicit function of the system states, as presented in Appendix 
A.2. 

The algebraic equations for the chosen output parameters were combined 
to obtain a matrix output equation in the form of Equation (4.23) using the 
output vector, 

y T ^ [ 7, n„ 0 r» #t 1 • (4-28) 

The numerical values of matrices C, E and F for each configuration appear in 
Appendix A.3. 

The pilot parameters, the time constants and D.C. gains were chosen to 
accurately describe the bandwidth limitation of the human pilot. A 
characteristic lag of 0.15 seconds was chosen to be consistent with other studies 
[8]. The resulting pilot filter equation used in the analysis is, 

*p =~ X P + j-1p (4-30) 

T p r p 

- where r p = 0.15 secs. 

The gust parameters were chosen to be consistent with a previous study 
using the B-l vehicle and the Dry den Gust Model (0], The gust equation used 
in the analysis is, 



(4.31) 



-9.4 

0.0 

+ 

1.0 

-0.0225 

-9.5 

0.0056 


V 


- for a g in radians and q g of unit intensity. 


Numerical Results 


The modal analysis method was applied to the eight configurations of the 
data base described in Chapter 3. The complete numerical results can be 
found in Appendix A.3. 

Consider the graphical results on Figures 4.2 - 4.9 which are the 
normalized relative magnitudes of the modal impulse residues for each mode of 
the vehicle due to pilot inputs. The normalization was done so that the residue 
magnitudes of the vehicle modes (not including pilot lag, i.e. phugoid, short- 
period and aeroelastic) sum to unity for each output. The equation used to 
accomplish this is Equation (4.32). 


R; 


i | norm 


Rj 

m 

mi 

j=« 


. (i = l m) 


(4.32) 


- where m is the number of vehicle modes. 


The absolute magnitudes for each mode can be obtained from the numerical 
results in Appendix A. 3. Since pitch atwtude and pitch attitude rate are 
outputs of primary concern in a pitch tracking task, the residue magnitudes 
associated with the rigid-body and total-elastic pitch attitude angles, (0 R and 
0 T ), and rates, (0 J{ and 0 T ), will be considered first. Clearly, 0 T and 0 T have 
more aeroelastic mode contribution than do and 0 R , which is as expected 
since 0 R is the rigid-body attitude angle and Oj is the total attitude angle 
including elastic deformation at the cockpit. 
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Now consider the results for 0 R and 0 r as the frequency of the first elastic 
mode is reduced as in configurations 1 through 3 (see Table 3.1 for reference). 
The residue magnitudes of the first aeroelastic mode (El) monotonically 
increases until, in Configuration 3, it is larger than the short-period modal 
residue! This indicates that, for Configuration 3, the rigid-body attitude 
response is dominated by the first aeroelastic mode! It is obvious that the use 
of a pure rigid-body analysis would be wrong and any model not including the 
effects of elastic modes would be inappropriate. 

The results also explain why Configurations 3 and 4, while having similar 
eigenvalue characteristics, have very different simulation result.', (see Table 
3.2). The aeroelastic modes in Configuration 4 do not dominate the attitude 
response (as they do in Configuration 3). The residue magnitude for the lowest 
frequency aeroelastic mode (E3 in this case) is not larger than the short-period 
residue. In other words, Configuration 4 has attitude dynamics which are 
dominated by a rigid-body mode and Configuration 3 has dynamics which are 
dominated by an aeroelastic mode. S.nce Configuration 4 acts more like a 
“rigid vehicle” than Configuration 3, the tracking performance for 
Configuration 4 is better than Configuration 3. However, the aeroelastic mode 
residue in Configuration 4 still contributes to some degradation in tracking 
performance. 

This approach can also be used to relate the rest of the tracking 
simulation results to the effects of the aeroelastic modes. The tracking errors 
(Figure 3.4) and the Cooper-Harper ratings (Figure 3.5) of the simulations 
agree especially well with the trends in the magnitudes of the impulse residues 
for total pitch attitude angle (0<]>). The configm tions with large tracking 
errors and poor pilot ratings have aeroelastic residue magnitudes which are 
larger than the rigid-body residue magnitudes in the 0? response. The converse 
is also true; the configurations with large aeroelastic residue magnitudes tend to 
have large tracking errors and poor pilot ratings. 

The graphical results can be used to bring attention to other aspects of the 
vehicle dynamics as well. Take, for instance, the plunge acceleration at the 
pilot station (n,). This parameter was, of course, of no importance in the 
fixed-based simulation, but would be of particular interest if the configurations 
were to be studied using a moving-base or in-flight simulator. The graphical 
results of the n, modal impulse residue magnitudes in Figures 4.2 - 4.0 indicate 
that ignoring aeroelastic affects when considering, for example, ride quality 
would be improper. The contribution of the aeroelastic modes is very 
important in the n, response of the vehicle for all configurations. 



Consider also, the flight path angle (q) response for Configuration 3. 
Ignoring aeroe astic affects in this case would give erroneous results since the 
aeroelastic residue magnitudes are significant compared to those of the rigid- 
body modes. 

Finally, the insignificance of the second and fourth aeroe’.jtic modes (E2 
and E4) in the pilot impulse response is clearly evident from 'he graphical 
results. The exclusion of these two modes in Yen’s (2) simulation study was 
tnerefore valid. 

The same type of trends i residue magnitude occur in the gust- 
disturbance impulse residue magnitudes, Figures 4.10 • 4.17. These results 
indicate that the aeroelastic modes contribute, in varying degrees, to the 
various vehicle responses due to an impulse input to the Dryden gust model, 
where an impulse input is the deterministic counterpart to “white" noise. Of 
particular interest are the results for rigid-body pitch-attitude-rate (0r). For 
Configurations 3, 4 and 6, the contribution of the aeroelastic inodes b very 
significant. One of the wing bendings modes, E2, has significant residue 
magnit udes compared to those of the rigid-body modes. Thb indicates that 
attitude tracking in turbulence would be similar for each cf these 
configurations in that the 0 R responses would be dominated by aeroelastic- 
structural modes. This implies that even though Configuration 4 had a 
satisfactory pitch attitude response in the simulation, added turbulence may 
result in significantly different and degraded performance. 

Tin modal analysis paints a different picture than the eigenvalue analysb 
presented in Chapter 3. Recalling the discussion in Chapter 2, one cun see that 
a* the frequencies of the structural modes are reduced, the interaction between 
the rigid-body modes and the aeroelastic modes increases. The resuM is that 
the residues associated with the structural modes and those associated with the 
rigid-body modes are modified and, as a result, alter the vehicle dynamics. If 
the residues of the aeroelastic modes become large enough to dominate the 
vehicle response, the aircraft no longer acts like a “rigid aircraft". In other 
words, the vehicle altitude response is not dominated by the characteristic 
short-period attitude dynamics. 

Since the residue magnitudes are a measure of *hc modal participation, the 
above argument indicates that when the impulse residue magnitudes associated 
with aeroelastic modes dominate those of the short-period inode, the vehicle 
performance degrades. The modal analysis results support this argument. 
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Gust Impulse Residue Magnitudes - Config. 6 
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Gust Impulse Residue Magnitudes - Config. 7 
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Ir. summary, the modal analysis method described iD this chapter has been 
used to attack the question of, “how and when do aeroelastic effects 
significantly afTect aircraft dynamics?” The analysis indicates that when the 
magnitudes of the modal impulse residues of the “aeroelastic modes” become 
the dominant residue magnitudes of the vehicle system for important outputs, 
the vehicle dynamics change significantly and may change in such a way as to 
result in “un-aircraft like” characteristics. In addition, the trends in the 
relative residue magnitude values for some outputs are closely related to the 
pilot ratings and tracking errors of the simulations. 

A drawback of using this modal analysis approach is that it is essentially 
open-loop in nature. Even though the modal analysis procedure considers some 
aspects of the pilot, specifically his limited bandwidth, it is still an open-loop 
analysis method. Since the pilot/vehicle system performance is really 
determined by the dynamics of the vehicle when the pilot closes the loop, a 
“closed-loop” or “pilot-in-thc-loop” analysis may give more insight into the 
effects ot the aeroelastic modes. The next chapter considers such an approach. 
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CHAPTER V 

CLOSED-LOOP ANALYSIS 


Complete flight vehicle system dynamics are dependent, not only on the 
aircraft dynamics, but also on the dynamics of the pilot and on how he 
interacts with the aircraft dynamics. Though the modal analysis method did 
consider the bandwidth limitations of the pilot, the method was still open-loop 
in nature. This chapter will apply a closed-loop analysis procedure to the 
configurations in the data base to study the effect of aeroelastic modes on 
closed-loop dynamics. 

The closed-loop analysis procedure that will be used here is an extension 
of the Neal-Smith procedure [10] which uses an optimal control model (OCM) 
of the pilot [11] in a pitch tracking task. This approach has, in the past, been 
applied to study the effect of flight control system dynamics on pitch tracking 
performance of fighter-type aircraft [12,13]. Since flight-control system 
dynamics and aeroelastic modes are both examples of higher order dynamics, 
there is reason to believe that this procedure may be useful in evaluating the 
effects of aeroelastic modes on the pitch tracking performance of the data base 
configurations. 

Before using this procedure to study the data base configurations, the 
procedure must be extended for application to flexible aircraft. This entails 
understanding the Neal-Smith methodology and applying the OCM to the 
Neal-Smith approach. This will be accomplished by briefly reviewing the work 
done by Neal and Smith [10] and by Bacon and Schmidt [12]. 


Neal-Smith/Bacon Methodology 

The investigation performed by Neal and Smith in the early 1970 ’s 
resulted in a criteria developed to expose problem areas in aircraft where the 
pilot was to perform a given task. Their criteria utilizes a closed-loop or 
“pilot-in-the-loop” analysis procedure. This pilot- in-the-loop approach was 
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used because of difficulties encountered in using existing open-loop handling 
qualities criteria and because of the truly closed-loop nature of piloted vehicles. 

The analysis method was based on the fact that the subjective pilot rating 
of a pitch- attitude task is primarily determined by how well the pilot can 
control the pitch attitude and how difficult it is to do so. Specifically, the 
analysis was performed using a compensatory tracking task model (i.e. the pilot 
only perceives the difference between the attitude of the aircraft and the 
commanded attitude), and by representing the pilot by a describing function 
consisting of a time delay and a lead-lag compensator (see Figure 5.1). The 
time delay accounts for perceptual delays and neuromuscular lags and the 
lead-lag compensation is used as a first order approximation of the pilot’s 
dynamic compensation in the task. 


Ik coiisidc.ing the closed-loop frequency response of the modeled 
pilot /aircraft system, Neal and Smith were able to relate the pilot's objective 
ratings to frequency response characteristics depicted in Figure 5.2. The 
resulting specifications also relate to the stated goals of actual pilots. 


For good performance, a pilot wants to be able to acquire the target 
quickly and predictably and with a minimum of overshoot and oscillation. 
“Quick acquisition of the target” implies that the pilot wants to achieve a high 
bandwidth. Neal and Smith also related minimizing overshoot to minimizing 


the closed-loop resonance peak, ( 


Or 


This inference conns from the 


relationship between system damping and the magnitude of the resonance peak 
in a second order system. By combining the two objectives, Neal and Smith 
concluded that, “pilot rating is a function of the compensation required to 
achieve good low frequency performance and the oscillatory tendencies that 
result.” (10.12,13) 


These objectives were related to the closed-loop analysis by defining the 
system bandwidth to be the frequency at which the closed-loop system phase 
lag reaches -1)0 deg.ees as illustrated in Figure 5.2. The pilot compensation 
was defined ns the phase of the resulting pilot describing function at the 
bandwidth frequency, excluding the efTect of the pure time delay, as shown in 
Equation (5.1). 


Ts(r,s+i) r s Aircraft 



Ncd-S-uith Model Structure 





Figure $.2 

Frequency Reeponee Specifications 
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PC ^ L 


J^BW T p, + [ 

jw BW T pa + I 


(5.1) 


- v/lmre PC *he pilot compensation. The closed-loop resonance peak is 
defined to be the maximum value of the magnitude of the closed-loop 
frequency response, (see Figure 5.2) 

The choice of parameters in the pilot describing function (K p , r, T pi , T pj ) 
were made to best satisfy a set of performance requirements. The requirements 
consist of - 


1) a specified time delay (r), 

2) a specified bandwidth characteristic of the task, 

3) a maximum allowable “low frequency droop", 

T pi 

-1) compensation (i.e. the value of — — ) required * .r a minimum value of 

* Pj 

resonance peak. 

Neal and Smith found that the parameters which satisfied these requirements 
resulted in a pilot phase compensation and closed-loop resonant peak that 
correlated with pilot rating. 

By plotting the value of resonance peak against PC (i.e. pilot 
compensation), Neal and Sni'th found regions in which aircraft with similar 
pilot rating were grouped. Figure 5.3 shows the Ncal-Smith result and the 
regions which correspond to the three levels of handlings qualities described in 
MIL-F-8785C (I). Level I includes aircraft having pilot ratings (Cooper- 
Harper' of 1.0 * 3.5, level 2 includes pilot ratings of 3.5 - 6.5 and level 3 ratings 
of 6.5 - 10.0. 

There are problems associated with the Neal-Smith method however. 
Th<»se problems lie in the difficulties associated with choosing appropriate 
frequency response specifications. For instance, it may be very difficult to 
determine the proper bandwidth frequency for an aircraft in a particular task 
without experimental data. Another problem lies in the choice of a maximum 
low frequency droop. Since the droop is only a relative measure of low 
frequency tracking performance, the choice of a maximum allowable value is 
rather arbitrary. Stiil fu?*h»r, the Neal-Smith method uses a compensatory 
task that was not representative of the actual task used in their flight tests. 
Bacon (13) extended the work of Neal and Smith to try to address these 
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problems. He applied an optimal control model of the pilot (OCM, [11]) to the 
closed-loop analysis method. The use of the OCM provides more flexibility in 
conducting the analysis by allowing the Neal-Smith criterion to be applied to 
other, more general, piloting tasks. It also eliminates the requirement of 
choosing the arbitrary frequency response specifications which are required to 
determine the pilot describing function. 

The optimal control model (OCM) of the pilot is based on the assumption 
that a well trained, highly motivated pilot choses his control input (u p ), subject 
to physiological limitations, in such a way that a quadratic cost function, 


J p = E 


lim t|; 
T— oo T 


/( 


X T Qx + ru* + g 


U p) 


dt 


( r »-2) 


is minimized. Here, Q and r are weightings chosen to reflect the task 
objectives and g is usually chosen to reflect physiological limits. 

Although details concerning the OCM can be found in [11], a brief 
description will be included here. Consider the block diagram of the OCM in 
Figure 5.4. The human perception characteristics that are modeled involve the 
pilot observations (y p ), passed through a pure time delay and contaminated by 
white noise of intensity V y . (see Equation 5.3) 


X P (t) = + liy(t-r) 


(5.3) 


X(t) = C p i(t) 

The solution to the stated optimal control problem yields a Kalman filter to 
estimate the delayed states and a least- mean-square predictor to obtain a 
current estimate of the states, ( £(t) ). The control law, obtained from 
minimization of the cost function J p , for a scalar u p , can be expressed as, 


= -K.i - u p 


(5.4) 



01 SUM AMCC S 
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- where K x is the optimal control gain matrix. The neuro-motor lag (r n ) 
Jesuits from including control rate (u p ) in the cost function J p and by 
weighting it so as to obtain a physically achievable value of r n , based tin man- 
machine experiments. 

As discussed in [12,13], when to OCM cost function is used to minimize 
tracking error, (0-0 c ), the resulting controller automatically minimizes low- 
frequency droop and resonance peak of the closed-loop system frequency 
response. This is an alternative to specifying maximum droop and determining 
the compensation required to minimize resonance peak in the Neal-Smith 
approach. In addition, the OCM will automatically determine the achievable 
bandwidth of the closed-loop system. Therefore, the Neal-Smith requirements 
are compatible with the OCM. The task for the analyst is now to properly 
apply t! e OCM. 

The proper application of the OCM involves, 

1) selecting a realistic pilot observation vector for the task (y p ), 

2) defining the cost function to be minimized (J p ) in the task, 

3) defining the command signal to be tracked (0 C ), 

•1) specifying the noise variances, observational thresholds and delays 

consistent with the human operator. 

By proper choice of pilot observations, cost function and command signal, the 
,anal> sis using the OCM closely reflects the inflight tracking task used by Neal 
and Smith. 

A subtlety discussed by Bacon [13] was associated with the almost- 
guaranteed stability of the OCM solution. Wi the cost function reflects 
minimizing tracking error, the resulting contro tries to make the closed-loop 
system act like a perfect tracker, that is, a system with a response-to-command 
transfer function of unity. As a result there is a trade-off between the low 
frequency droop and resonance peak. That is, an aircraft that would actually 
lead to oscillatory tendencies and a significant resonance peak in the Neal- 
Smith analysis, would yield an OCM solution that would sacrifice low 
frequency performance for stability. This is analogous to the pilot being less 
aggressive and tracking the target so that the oscillations would not occur. 

This piloting strategy tends not to expose the tendency of the aircraft to 
oscillate. 

Bacon argued that, since oscillation occurs from “suboptimal” piloting 
strategy, the OCM controller would do a better job of tracking than a real 
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pilot. He further argued that by increasing the forward path gain, one could 
approximate an aggressive pilot. An aggressive pilot would try to obtain better 
low frequency performance at the expense of high frequency oscillations. This 
type of piloting strategy would therefore expose the oscillatory nature of an 
aircraft. Hence, in Bacon’s approach, the forward path gain (i.e. K p in Figure 
5.5) was adjusted so that each vehicle configuration led to a maximum low 
frequency droop, similar to the Neal-Smith method. The adjustment exposed 
the oscillatory aircraft by increasing the resonance peak of such aircraft. 

Figure 5.5 illustrates this effect. 

By plotting the gain- adjusted resonance peak against the pilot phase 
compensation (obtained from the OCM) at the bandwidth frequency (as 
illustrated if Figure 5.6), Bacon obtained a plot analogous to the one obtained 
by Neal and Smith. Figure 5.7 presents Bacon’s results for the Neal-Smith 
Configurations which can be compared with Neal-Smith ’s original results shown 
in Figure 5.3. Just as in the Neal-Smith criterion, aircraft with similar pilot 
ratings group together on Bacon’s plot. 

The choice of proper bandwidth is not necessary in the Bacon method and 
is replaced by choosing a weighting in the cost function which results in a 
reasonable neuro-motor lag (r n ), which is a natural physiological limit. Bacon’s 
[13] results also indicate that the closed-loop system bandwidth, a result from 
the OCM analysis, correlates with subjective pilot rating. In fact, this 
relationship has been suggested as an additional criterion for measuring the 
quality of the vehicle dynamics.[l2,13] Figure 5.8 illustrates the relationship 
between closed-loop bandwidth and pilot rating for the Neal-Smith 
Configurations. 

A disadvantage of both methods is the need to choose au arbitrary 
maximum low frequency droop. As an extension of the Neal-Smith/Bacon 
method, an alternate way of consider*!.. g oscillatory tendencies will be 
presented. This new variation of the the Neal-Smith/Bacon method will be 
used to consider how aeroelastic modes affect the closed-loop characteristics of 
aircraft. 


Extension Of Neal-Smith/Bacon Methodology 

Id an attempt to make the analysis procedure independent of an arbitrary 
choice of the maximum low-frequency droop, an alternate method is proposed. 
Pilot induced oscillations (PIO’s) usually occur with aggressive pilot behavior. 
If the pilot "backs-off” (i.e. reduces his aggressiveness), the oscillations 
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disappear. A poor aircraft may have characteristics which produce PIO’s, with 
just slight increases in aggressiveness. It is this characteristic which Bacon’s 
gain adjustment exposes. 

A first order approximation to pilot aggressiveness is the DC gain of the 
pilot describing function. A slight increase in this gain from the OC'M 
approximates an increase in pilot aggressiveness. If the increase in closed-loop 
resonance peak which results from the increased pilot gain is relatively large, 
one could conclude that the aircraft is sensitive to pilot aggressiveness or, in 
other words, has oscillatory tendencies. 

This argument implies that a gam sensitivity procedure can be used to 
expose aircraft with oscillatory tendencies. Since the OCM optimizes the 
controller design in such a way that low frequency droop is sacrificed for 
resonance peak, excessive droop results for configurations with oscillatory 
tendencies. Using these ideas, a gain sensitivity parameter is defined to be, 


SP ^ DROOP x 



(dB). 


(5.5) 


The term. DROOP, is the magnitude of the “droop” for the case in question, 
obtained directly from the closed-loop OCM analysis (see Figure 5.2). The 
other term on the right hand side is the relative gain sensitivity which is 
determined by calculating the change in resonance peaks obtained using the 
basic pilot (model) gain and that obtained using a perturbed gain, £K. The 


sensitivity is the ratio of the change in resonance peaks, 
gain difference (AK). * 


ikl* 


the 


To justify the validity of using SP as a measure of oscillatory tendency, 
consider Table 5.1. This table presents the resonance peaks of the 
configurations from the Ncal-Smith study and the Bacon study and the values 
of the SP calculated for the same configurations. 


t The AK that was used to obtain the numerical results was determined by perturbing 
the pilot DC gain by approximately ten percent. 


W 
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Table 5.1 

Comparison of Resonance Peak and SP Values 
Resonance Peak (dB) 


Config. 

Neal-Smith 

Bacon 

SP (dB) 

mm 

7.0 

7.10 

133 

19 

0.5 

1.86 

0.40 

mm 

2.0 

4.84 

1.13 

ID 

00 

1 83 

0.30 

IE 

00 

3.50 

0.73 

2A 

3.0 

4.07 

1.50 

2B 

7.0 

11.37 

2.44 

2C 

2.0 

1.20 

0.03 

2D 

2.0 

1.24 

0.84 

2E 

3.5 

3.28 

1.40 

2F 

2.5 

3.00 

1.20 

2G 

6.0 

0.25 


2H 

3.0 

2.50 

0.81 

21 

7.0 

6.36 

160 

3A 

-1.0 

0.68 

0.28 

4A 

10.0 

10.17 

2.26 

5A 

>12 

18.21 

3.71 

6C 

1.5 

1.25 

0.35 

7C 

00 

0.77 


8A 

0.0 

0.646 

WSfiM 
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The trends between the three parameters, for the most part, agree well with 
each other. Furthermore, Figure 5.9, when compared with Figure 5.7. indicates 
that SF is an analogous measure of oscillatory tendencies and can therefore be 
used instead of resonance peak in the analysis. 


Application Of The Neal-Smith/Bacon Analysis To 
The Data Base Configurations 

In order to apply the Neal-Smith/Bacon analysis method, a clear 
understanding c * the similarities and differences between the type of 
configurations studied by Bacon and the flexible aircraft of the data base of 
Chapter 3 is necessary. The aircraft used in Bacon’s study were some of the 
configurations used in the Neal-Smitb study, and represent basic airframe 
dynamics with control system dynamics added. The basic aircraft dynamics 
that were analyzed included only the short-period me de. The added high order 
modes representing control system dynamics included a real pole, a real zero 
and a second-order, oscillatory mode. Figure 5.10 shows the basic airframe 
plus flight control system (FCS) dynamics in the pitcb-attitude-rate-to-stick- 
deflection transfer function. The short period dynamics are determined by Tj a 
and by u>, p and ( Ip and the r CS dynamics are dt’.ermined by r,, r 2 , u 3 and f 3 . 

The flexible aircraft of this study also have higher order modes but they 
correspond to aeroelastic effects and not FCS effects. The dynamics of the 
flexible aircraft have already been discussed in Chapters 2 and 4 and it will 
suffice to summarize them with Figure 5.11. Here the rigid-body dynamic 
parameter: are T #i , T #a , u> ph , f ph , u> tp and f fp , and the aeroelastic effects lead to 

<Tj , £ (i = l m) and u>j, ^ (i = l,...,m) where m equals the number of structural 

modes included in the vehicle model. 

An important step in the analysis is to decide on an appropriate cost 
function (7 p ). In the case of flexible aircraft the pilot sees or senses total 
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Figure 5.II 

of Ftatible Configuration* 
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tracking error (c T ~ 0j- 9 C ), with structural effects included. Choosing total 
error as the minimized parameter, however, may not correctly reflect the pilot’s 
objectives. Tim pilot comments from the simulation of the data base 
conugurations suggest that the pilot tried to track rigid-body error. 

The following quotes are typical of the pilot comments that resulted from 
the simulation study (2). 

For Configuration 8, the pilot comments included - 

“ more oscillation involved due to elasticity apparently, but it was high 
enough frequency that it was easy to ignore that and simply to fly the 
rigid body ...” (2) 

For Configuration 7 the pilot comments included - 

“ it was not too difficult to ignore (the oscillation) and to fly the rigid 
(body) ...” [2] 

These comments indicate that the pilot places more emphasis on keeping the 
rigid tracking error low than on minimizing total (displayed) tracking error. 
(Also see (12,14).) Therefore, the appropriate cost function for the flexible 
aircraft is, 


•W 


= E 


lim 

T-ooT 




(5.6) 


- where — t? c ) and g is chosen to obtain the desired r n . 

Bacon has shown (12,13) that the choice of r n takes the place of bandwidth 
in the Neal-Smith method. The value of r n is chosen to reflect pilot 
aggressiveness in the tracking task and determines the bandwidth of the 
closed-loop system. As r n increases, by increasing g, the bandwidth decreases. 
Low r n , which represents ‘ aggressive behavior”, results in a fast closed-loop 
system and so a high bandwidth. Therefore, to obtain the maximum possible 
•'andwidth, r n should be set at the lowest physically possible value, which is 
usually coi. idered to be 0.1 seconds. (8) The value of g used in this study to 
obtain a r n of approximately 0.1 seconds was, 
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g = 0.0040 . 


The analysis includes the v hide dynamics, the pilot observations and the 
command signal dynamics. These factors must be chosen to be consistent with 
the task. The complete pilot observation vector therefore includes 
e T and < T , 0 T and ff T . These four parameters are naturally displayed to the 
pilot in the tracking task. 


The intermediate output of the analysis consists of the controller gains, 
closed-loop eigenvalues, rms values of the inputs, states and output parameters, 
cost function values, optimal estimator gains and, most importantly for this 
application, frequency responses for selected transfer functions. By combining 
the transfer functions properly, the desired closed-loop transfer function 
frequency response, similar to that of Neal and Smith, can be formed. 


Consider the block diagram of the tracking task in Figure 5.12. The 

0r(s) 

closed-loop transfer function of interest in this study is — ■ , since rigid-body 

0c( s ) 

attitude, 0 R , is what the pilot “cares about” in rating the performance of the 
aircrdft. This transfer function can be written as, 


W = H(s)G|(s) 

*cM 1 + H(s)G,(s)G 2 (s) 


e T (s) 


1 + 


0t( 8 ) 

<t( s ) 


(5.7) 


0 R (s) 0 c (s) 

The transfer functions of interest are therefore — — and — —• Table 5.2 

€ X («) < T (s) 

summarizes the numerical values used to obtain the desired results from the 
analysis. 


The vehicle dynamics (Equations 4.21-4.23) are the same as those used in 
the open-loop nrodal analysis. These consist of the A v and B matrices of the 
data base configurations. The D matrix is zero however since gust disturbance 
dynamics is not considered in this analysis. 
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Table 5.2 

Summary of Closed-Loop Analysis Inputs 


Observational 

Thresholds 

0R • 0Ji 
0R > 0J t <T 

8.7 x 10~ 4 rad * 

3.1 x 10' 3 — * 
sec 

Variance of 

f T , «T t 

-20 dB * 

Sensor Noise 

% - 

-6dB 

Attention Allocation 

f T , «T' 

0.245 * 

(V A.A. i = 1-0 ) 
i 

0R> ^R 

0.01 


* consistent with previous work [13] 


The remaining requirement, is the command signal dynamics. A command 
signal (#(<), which accurately represents a challenging pitch tracking task and 
approximates the tracking task used by Neal and Smith, is defined by Equation 

5.8. 


0 C + 0.5 <) c + O.250 c = w(t) (5.8) 

Here. 0 C is the commanded attitude and w(t) is zero mean (iaussian white noise 
of intensity Y w . The intensity of the white noise was chosen to result in an 
rtns value for Oq- of approximately three (3) degrees. 

The resulting, model-compatible, state variable representation has the 
form, 




A r 

0 


0 1 

E e 

• 

X.OCM ~ 

0 

Ay 

'S-OCM + 

ByP 

Of 

• 

where. 







T — 

X-OCM “ 



1 J | 




(5.0) 


(5.10) 
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Here, A c and E c are the matrices resulting from the state variable 
representation of the command signal. 


Closed-loop evaluation of the model yields the desired frequency responses 
and the pilot describing function frequency response. The desit ed closed-loop 

0 R (s) 

frequi response, namely - is obtained by manipulating the frequency 

0 c (s) 

responses according to Equation 5.7 at selected frequencies. That is, 


0 c (j^) 


w 

M s ) 

^ 0t( s ) 
c x (s) 


S ~ jut 


(5.11) 


The frequency responses that result from the closed-loop analysis of the eight 
data base configurations, can be found in Appendix A.5. An example of the 
frequency responses is shown in Figure 5.13. The “Purdue Pilot” frequency 


response corresponds to H(s) or 
response corresponds to Gj(s) or 


M 

£ X (s) 

*r(s) 


and the “Aircraft (O.L.)” frequency 


as depicted in Figure 5.12. The 


“Aircraft Plus Pilot (O.L.)” frequency response corresponds to H(s)G t (s) or 
0 R (s) 

— — and "Aircraft Plus Pilot (C.L.)” corresponds to the frequency response 
c x (s) 

0r( s ) 


for 




Numerical Results 

The closed-loop system frequency response properties; bandwidth, droop, 
resonance peak and sensitivity parameter; pilot phase compensation at the 
bandwidth frequency and pilot rating are summarized in Table 5.3 for each of 
the data base configurations. 

First examine the trends in pilot rating with closed-loop bandwidth. 
Figure 5.14 is a plot of pilot rating (PR) versus bandwidth frequency (u> BW ) for 
the eight data base configurations. Though the number of data points is 
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Table 5.3 

Summary of Closed-I top \nalysis of Data Base Configurations 


Case 

rad 

^BW 

sec 

DROOP df 

1 

1 

max dB 

SP dB 

PC de* 

PR 

1 

2.29 

i.O 

msm 

M 

-58.9 

1.6 


1 80 

1.51 

HI 

ESI 

•66.6 

2.0 


0.07 

8.45 

6.99 

7.59 

-96.02 

5.9 


1.77 

1.93 

-0.11 

0.78 

-16.08 

3.1 


1.53 

2.14 

0.10 

0.49 

-68.9 

2.0 

6 

0.10 

6.94 

9.59 

wm 

-72.3 

6.7 


1.70 

1.56 

1.02 

rm 

-60.55 

2.3 


2.00 

1 23 

1.69 

m 

-53.3 

1.9 


limited, the trend is consistent with that of (12,13). The configurations with 
low closed-loop bandwidth have poor pilot ratings and the configurations with 
relatively high bandwidth have better pilot ratings. 

Next, consider the closed-loop system parameters of Neal and Smith. 
Figure 5.15 is a plot of “Noal-Smith like” criteria for the eight data base 
configurations. However, in place of resonance peak, the sensitivity parameter 
(SI*) is used as a measure of oscillatory tendencies. Note that the 
configurations with similar pilot ratings (Pit) are grouped together and the 
configurations with poorer ratings (i.e. configurations 3,4 and 6) are distinctly 
separated from the better aircraft. Also, two of the configurations (i.e. 3 and 6) 
have relatively large values of SI*, indicating oscillatory tendencies. This 
oscillatory nature is also noted in the pilot comments from the simulations (see 
Table 3.2). Notice that the value of SP for Configuration 4 indicates that its 
poorer performance is not due to oscillatory tendencies. The pilot 
comp 'iisation (PC), though, indicates that the pilot has to supply more lead for 
the best performance, or in other words, the aircraft response is sluggish. This 
sluggish nature is also noted in the pilot comments from the simulations (see 
Table 3.2). 

Therefore, the analysis not only grouped aircraft with similar pilot ratings, 
but it also exposed response characteristics that contribute to degraded 
performance. Though there is not enough data to determine boundaries 
defining the three handling qualities levels, the trends tend to imply their 
existence. The implication is that this closed-bmp analysis might be abic to 
identify when ae.*oflailic effects significantly affect the dynamics of flexible 
aircraft. That is *he .'.'cal-Smith/Ha* - a analysis, appropriately utilized, may 
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appl; to large flexible aircraft as well as small aircraft with added control 
system dynamics. More specifically, the results indicate that the data base 
configurations with poor tracking performance received the poor results because 
of sensitivity to forward path gain (used to approximate pilot aggressiveness) 
and indicates oscillatory tendencies. 



CHAPTER VI 
CONCLUSIONS 


The objective of this study was to investigate when and how structural 
effects (especially dynamic aeroelastic effects) affect the dynamics of aircraft. 
Two analysis methods, an open-loop modal technique and a pilot-in-the-loop 
method, were used to see how aeroelastic modes affect the dynamics of aircraft 
in the longitudinal axis. Both procedures were applied to a family of aircraft 
which exhibit considerable aeroelastic effects. 

The results of the modal analysis indicate that when the magnitudes of 
the modal impulse residues of the aeroelastic modes become large compared to 
the residue magnitudes of the rigid-body modes for important outputs, the 
dynamics can change significantly and in such a way that the handling qualities 
of the vehicle may be degraded. In addition, the trends in impulse residue 
magnitudes for some inputs are closely related to the trends in pilot ratings of 
the configurations from the fixed based simulation. 

The pilot-in-the-loop analysis verifies that as the frequencies of the 
aeroelastic modes decrease, the performance of the vehicle tends to degrade. 
More specifically; as the structural vibration frequencies were decreased, the 
sensitivity of the closed-loop system to perturbations in forward path gain 
increased. This effect was demonstrated by plotting feed-forward gain 
sensitivity (SP) versus pilot compensation (PC) in a tracking task. It was also 
shown that the bandwidth of the closed-loop system correlates with the 
subjective pilot ratings and those configurations with lower structural 
frequencies tend to have lower closed-loop bandwidths. These results indicate 
that reduced aeroelastic mode frequencies can cause degraded handling qualities 
which may appear in the form of oscillatory tendencies and sluggish response. 

In conclusion, dynamic aeroelastic effects can definitely contribute to 
degraded performance of aircraft in the longitudinal axis. The aeroelastic 
modes contribute to poor performance primarily by, 1) introducing dynamic 
effects of their own in the form of high frequency oscillations, and 2) modal 
interaction which alters the dynamics cf the rigid-body modes. In addition, 



these effects can occur when the aeroelastic mode frequencies are still several 
times higher than the frequencies of the rigid- body modes! As a consequence of 
these effects, aeroelastic modes should be taken into account for vehicles that 
exhibit these dynamic aeroelastic effects. 

Future work in this area should include expanding the data base. With a 
larger set of configurations to study, the analysis methods developed here can 
be applied to obtain more conclusive results which may lead to quantitative 
rules for specifying handling qualities for flexible aircraft. For example, it may 
be possible to define handling qualities boundaries in the SP versus PC (i.e. 
sensitivity parameter versus pilot compensation) plot from the pilot-in- 
the-loop analysis. The boundaries would divide the plot into three regions 
which correspond to the three handling qualities levels. Also, the analysis 
methods developed here should be extended to study lateral-directional 
dynamics in order to understand the problem more completely. Finally, since 
it has been shown that aeroelastic modes can be important, future work «hould 
be aimed at developing control synthesis techniques that utilize the modal 
techniques, either directly or indirectly, to gain insight into the consequences of 
aeroelastic effects. Such techniques might address restoring excellent handling 
qualities to vehicles with poor handling due to dynamic aeroelastic effects. 
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Appendix A.1 

Scaling Transformation for Mode Identification 


The aircraft states are scaled so that the elements of the eigenvectors have 
comparable units. This is done so that the eigenvectors can be used to simplify 
the task of identifying the modes of the system. That is, aid in determining 
which eigenvalues are associated with, for example, the short-period mode or 
one of the aeroelastic modes. 

The scaling of the system states is accomplished by means of a similarity 
transformation applied to the vehicle state variable model of the form, 

i = Ay + Bu 


(All) 


y = Cjl + Du- 

Consider a flexible aircraft in the longitudinal axis. The following state vector 
definition is representative of such an aircraft. 

y T ^ [ u, a, 9, if, ri ] (A. 1.2) 

In the longitudinal axis, pitch angle and pitch rate are two pertinent 
dimensions. The vehicle states can be scaled so that all of them are 
nondimensional or can be physically interpreted as angles and angular rates, 

(i.e. units of radians and radians per second). The forward velocity 
perturbation, u, can be divided by the cruise velocity, U 0 . The generalized 
elastic deflection, tj, can be multiplied by the mode slope, 0' , which makes the 
the state physically analogous to elastic pitch angle with units of radians. This 
is evident when considering the equation for total-elastic pitch angle, 
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i=l 


(A. 1.3) 


Similarly, the generalized rate, i), can be multiplied by the mode slope. The 
result is that the state becomes analogous to elastic pitch rate with units of 
radians per second. The rigid-body pitch attitude, 9, pitch rate, 0, and angle of 
attack, a, are expressed in radians and so do not need to be scaled. 

For the model and the scaling factors described above, the similarity 
transformation can be defined to be, 


“■ 0000 0 
Uo 


T = 


0 1 0 0 0 0 

0 0 1 0 0 0 

0 0 0 1 0 0 

0 0 0 0 0 ' 0 

0 0 0 0 0 ^' 


The transformed state vector is defined by, 


£ = T*. 


(A. 1.4) 


(A. 1.5) 


Applying the transformation to the vehicle model in Equation (A.1.1) results in 
the scaled system, 

i = TAT -1 £ + TBu, 


(A. 1.6) 


X — CT *£ + Du- 


An important property of a similarity transformation is that it has no 
affect on the eigenvalues or residues of the original system. This property 
allows the scaling transformation to be applied to the vehicle model without 
altering the results of the modal analysis. 



Therefore, by applying the scaling transformation described above, the 
units of the eigenvectors can be adjusted to make identifying the modes of the 
system easier. In addition, this can be done without affecting the results of the 
modal analysis procedure. 
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Appendix A.2 

d,ui Function of the Vehicle States 


Consider the longitudinal state variable model of an aircraft, Equation 
(A.2.2), with the following state variable definition. 


2 t. T i I u, a, 9, 9,r), t)} 


(A.2. 1) 


i = Ax. + Bu 


X = Ci + Du 


(A.2.2) 


The plunge acceleration of an aircraft (n f ) is described by the following 
expression, 


n,(t) = - 

g 


m 


I'oim + i.«n) - v; ah, I 'MU 


i = ! 


.(«■=) 


- where g = gravitational acceleration, (ft/sec 2 * 

U 0 = cruise velocity, (ft/sec) 

1, = distance from c.g. to cockpit, (ft) 

<b x = mode shape of i ,h acroelastic mode, (ft) 
m - number of acroelastic modes 


(A. 2.3) 


The other parameters in Equation (A.2.3) can be expressed in terms of the 
states of the aircraft. 
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The flight path angle is defined as, 


1<t) A*t)-o(t). 


(A.. :) 


Therefore, 

^t) = 0(t)- a(t). (A.2.5) 

Note that 0(t) is a state of the vehicle but a(t) is the time derivative of the 
vehicle state a(t). Note also that the derivative of the angle of attack with 
respect to time can be written as, 


o(t) = A^x. + B^u (A.2.6) 

- where A^ and B q are the rows of the matrices A and B, respectively, 
associated with the scalar equation for o(t). 

Similarly, 0(t) and i/(t) are the time derivatives of the states 0(t) and i)(t). 
Therefore. 


0(t) = Ajg + Bju 

(A.2.7) 

»?( t) = Ajjg. + B^u 

(A.2.8) 


- where Aj and Bj, and A- and B-^ are the rows of A and B associated 
with the scalar equations for 0(t) and »j(t), respectively. 

Using Equations (A.2.S), (A.2.6), (A.2.7) and (A.2.8), the expression for the 
plunge acceleration can be written as, 

Uo(«M- Ail + B iU ) + I. U,1 + r.j») 




s 


i 
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Note also that 0(t) can be written as, 


(A. 2.0) 


0(t) = A** + BjU (A.2.10) 

• where A t and B^ are the rows of A and B associated with the scalar 
equation for 0(t). 

Therefore, 





UoJa^-A^Jx. + u 0 (B r B o )u 


+ I, (a>jl + B,u)- _E^(a^ + B-u)J 


(A.2.11) 


By grouping the terms multiplying x. and u. simple expressions for the rows of 
C and D associated with the scalar n, output equation can be formed. 


C 


n. 


I 

s 


D 


n. 



L 'ok-A,j + i, a, - £ ^a, 

' 9 i= I 

u.(Ur®i) + 4%-£*V 


(A.2.12) 


(A.2.13) 


• wheie C n< and D n< are the row of C and D associated with the scalar n t 
equation. 

This method of determining the proper coefficients for the C and D 
matrices associated with n, can be implemented numerically very easily by 
using a transformation row vector, X. The definition of the transformation 
depends on the state vector for the system. For the state definition in 
Equation (A.2.1), the transformation vector has the following form, 
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X 4(0, 


Uo Uo 
g ’ g 



(A.2.14) 


Thus, 


t 


C n . = X A 

and, 

D».=XB 


(A.2.1S) 


V 


(A.2.16) ! 
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